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The representative measurements of the second harmonic in ac magnetization complemented by 
neutron depolarization have been performed for single-crystal Lao.ssSro.isCoOs in the temperature 
range 97 K< T < 230 K, where occurrence of a small fraction (< 10“®) of nanoscale ferromagnetic 
clusters (FMC) has been found. Magnetic, geometrical and dynamical parameters of the FMC 
system have been evaluated in the temperature range T < 140 K, where superparamagnetic regime 
installs, by means of the formalism involving the Fokker-Planck equation (FPE). With lowering 
the temperature, the amount of clusters fraction, the cluster size and magnetic moment along with 
its diffusion relaxation time strongly increase, each in its own temperature interval. Below 130 
K, FMC contribute essentially to the total linear magnetic susceptibility. The damping factor of 
the order 10“^ proves the importance of precession in thermal relaxation of the cluster magnetic 
moment. The FMC are a precursor of long-range ferromagnetic correlations seen below 100 K 
with neutron-scattering techniques. The employed technique supplemented with FPE-based data- 
treatment formalism is a novel method for studying superparamagnetic systems. 
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I. INTRODUCTION 


Hole-doped lanthanum cobaltites are known to exhibit the tendency to magnetoelectronic phase separation (MBPS) 
in the wide range of doping and temperature M- The parent compound LaCoOa in the ground state is the 
diamagnetic insulator. The 3d orbital of Co^+ ions is split by nearly cubic crystal field, slightly exceeding intra-ion 
exchange, in such a way that t 2 g orbital triplet is completely filled with six electrons, forming the low-spin (LS) state 
with the spin S' = 0, whereas the empty Cg orbital doublet lies 10-12 meV higher. At elevated temperatures, LaCoOa 
becomes paramagnetic due to the thermally induced transition to higher spin states. 

The hole doping with alkali-earth ions Sr^+ gives rise to nanoscale MBPS in the form of spin polarons and/or 
larger spin clusters depending on temperature and doping, as established from neutron diffraction, nuclear magnetic 
resonance (NMR), Mossbauer spectroscopy Q, small-angle neutron scattering, inelastic neutron scattering, x-ray 
spectroscopies, heat capacity, mag netometry, ac magnetic susceptibility (MS), thermopower, transport and m^neto- 
transport measurements (Refs. andQ and a wide bibliography therein). On the x-T phase diagrams P andQ, the 
MBPS state covers a broad area below the paramagnetic phase. The transition temperature increases with doping, 
so that, e.g. at x = 0.15, the spin-polaronic state extends up to T si 250 K [T|. At lower temperatures somewhat 
below 60 K, a highly inhomogeneous state develops P and[l], its particular form depending on x. Bor the doping less 
than the characteristic value Xc = 0.17, pronounced hole-rich and hole-poor regions coexist in the form of insulating 
spin glass, whereas in the doping range Xc < x < 0.25, the spin-polaronic state evolves upon cooling followed by 
the long-ran ge- order ferromagnetic (PM) region eventually passing to the metallic phase of strongly inhomogeneous 
ferromagnet P . According to other sources (Ref. Q and references therein), the insulator-metal transition occurs in 
the interval 0.18-0.22, with the percolation onset at Xc = 0.18. 

The nature and temperature evolution of MBPS close to Xc has been of keen interest. In the insulator-state vicinity 
of Xc, large-scale magnetic inhomogeneities in single crystals have been known to occur at temperatures not exceeding 
100 K implying spin-polaronic MBPS at higher temperatures P. An extensive neutron scattering study of the single¬ 
crystal Lao.ssSro.isCoOa revealed coexistence of the regions with long-range and short-range order at temperatures 
T < 100 K p. The former were identified as PM clusters of the diameter larger than 7 nm and the latter as referring 
to smaller, less than 2 nm, entities. Under external magnetic field, large clusters formed at the expense of small ones. 
After switching off the held, some part of the large clusters still persisted. The small-angle neutron scattering assisted 
by magnetotransport measurements for the same composition also showed spontaneous MBPS developing below 150 
K in the form of small PM clusters embedded in a non-PM matrix reaching the size 1.5-2.5 nm at low temperatures 


It is generally accepted that, similarly to manganites, in doped cobaltites double exchange interaction between 
Co^’*' and leads to formation of PM clusters with metallic conductivity with higher density of charge carriers 

than in the surrounding matrix. Increase of the Coulomb energy with the cluster growth is compensated by kinetic 
energy of the carriers up to some critical cluster size, above which the clusters become unstable |1(TI - [18I| . 

The MBPS occurs in manganites, as well, revealing a certain resemblance to that in cobaltites at a close level of 
doping, in compounds both conserving their insulator state and exhibiting insulator-metal transition below the Curie 
temperature [Ij andll^. In manganites, however, the increase of paramagnetic fluctuations upon cooling resulting 
in the PM phase transition impedes and restricts in temperature observation of PM clusters [ij]. The latter is easier 
for cobaltites as their matrix is only slightly paramagnetic down to lower temperatures. 

In polycrystalline LaSr cobaltites, PM clusters occur in the wider doping and temperature ran ges, as crystallite 
boundaries favor their formation. Prom linear and non-linear magnetic susceptibility measurements |I6l | , PM clusters 
were found to exist at the dopings x > 0.25 and to appreciably contribute well above 100 K even in the insulating 
state for x = 0.1 and 0.15. 

Previously, the non-linear magnetic response was found to occur in single-crystal Lao.ssSro.isCoOa in the wide 
temperature range below T Ri 213 K [la. The preliminary qualitative analysis has shown that, below this temperature, 
a small fraction of PM clusters emerges, exhibiting nonsuperparamagnetic behavior characterized by considerable field 
hysteresis. Localization of the clusters was attributed to preferable sites formed by slight doping heterogeneity and 
local oxygen nonstoichiometry. The concentration of clusters increased weakly with lowering the temperature. On 
cooling below T Ri 147 K, homogeneous nucleation of PM clusters developed exhibiting superparamagnetic (SPM) 
behavior with only a small field hysteresis. Starting from 135 K, the concentration of SPM clusters increased intensively 
in a narrow temperature interval. 

However, the former studies [l3, [H, and provide quite limited, indirect, and approximate information on PM 
clusters being based on the analysis of the response signal itself, viz., extrema positions, signal amplitudes, field 
asymptotes, etc. In this report, we specify the previous observations and present a detailed complete quantitative 
characterization of the cluster system involving a thorough mathematical formalism for data treatment especially 
focusing on the temperature region below 140 K. Bvolution of clusters was traced down to 97 K at magnetic fields not 
exceeding 300 Oe. From representative data on the second harmonic of longitudinal magnetic response complemented 
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by neutron depolarization results, a large set of parameters characterizing geometrical, magnetic, and dynamical 
features of the system of FM clusters in SPM regime was determined using the formalism based on the solution of the 
Fokker-Planck equation (FPE) [20l - [^ . With these parameters, we visualize the above mentioned successive phases 
of the cluster MBPS, specify the temperature boundaries of the MBPS stages, and outline two different regimes of 
homogeneous nucleation alternating each other at 115 K. 

The observed large SPM clusters are suggested to be (i) a precursor of the FM state occurring below the insulator- 
metal transition at the doping x > Xc and (ii) an onset of large FM clusters revealed by Phelan and collaborators 
with neutron scattering techniques at lower temperatures, T < 100 K 

The technique employed in this research, though not being conventional, proved its high efficiency in studying large, 
of the order 10 nm, FM particles [13) EM andEM|. Such objects appear to give the main contribution to generation 
of the second harmonic at magnetic fields of the order 10-100 Oe due to strong nonlinearity of their response in weak 
fields. In this study, the signal from the SPM fraction occupying the sample volume > 10“"^ was reliably measured 
and the parameters characterizing the SPM system were obtained with sufficient accuracy. The ac-field frequency 
matched the most informative frequency range where the real and imaginary parts of the second-order response signal 
were comparable. 

In Sec. 11, relevant experimental details both for ac and neutron depolarization measurements are presented. Section 
III acquaints with the principles of extracting the quantities characterizing the cluster system and explains how some 
more parameters can be evaluated involving additionally the neutron depolarization data. In Sec. IV, the neutron 
depolarization data and the non-linear response for some typical temperatures are presented and shortly described. 
Temperature dependencies of the characteristic parameters are analyzed and a scenario of the cluster-system evolution 
is discussed with the separate emphasis on the SPM dynamics. The linear ac susceptibility retrieved from the non¬ 
linear data is compared to the measured linear response. Section V is devoted to quite a special problem concerning 
the grounds of SPM dynamics. Experimental data on the non-linear ac response appear to be representative enough 
to distinguish between Landau-Lifshitz and Gilbert ansatz underlying the Fokker-Planck formalism and to make a 
choice in favor of the latter. In the conclusion section, the main results are summarized. The appendix contains the 
formalism for ac nonlinear data treatment. 


II. EXPERIMENTAL DETAILS 

The manufacturing and certifying of the specimen are described in detail in Ref. fl3 | . The feed rod for the sample 
was fabricated with the standard routine based on solid-state synthesis. The single-crystal Lao.ssSro.isCoOa was 
grown by the floating-zone technique using radiative heating under an oxygen pressure of 1 bar with a typical growth 
rate of 1 mm/h. The crystal was found to be single-phase, pseudocubic with the slight rhombohedral distortion at 
room temperature (the space group i?3c), as confirmed by x-ray diffraction. The element content was monitored by 
an x-ray microanalyzer with the relative accuracy 3%. 

The second harmonic of magnetization M 2 was measured in p arallel dc and ac magnetic fields on the setup 
approved in a good deal of magnetic measurements (Refs. EMEMEMlHi andl^ and references therein). The dc field 
H was scanned symmetrically with respect to zero field within ±300 Oe with the round-up cycle 0.14 s. A frequency 
and an amplitude of the excitation field were uj/2tt = 15.65 MHz and h = 0.8 — 14 Oe, respectively. The latter 
was optimized to obey the condition M 2 oc h^, necessary for the preliminary comparative data analysis, although 
not mandatory in the FPE data-treatment formalism. Both phase components of the signal, ReM 2 and ImM 2 , were 
simultaneously recorded as functions of the dc field in the temperature region 97 K < T < 230 K. The sample 
temperature was stabilized with an accuracy of 0.2 K. 

The neutron depolarization techni que is known to be a powerful mesoscopic technique to study ferromagnetic 
clusters in a phase separated system |28l - l3l| . The neutron magnetic moment interacts with the internal magnetic 
field of an FM cluster. As the internal fields of different clusters are thermally disordered, the neutron beam, passed 
through the cluster ensemble, depolarizes. The contribution of spin polarons to neutron depolarization via forward 
scattering is negligibly small, due to their small size. The scattering on structural inhomogeneities also does not affect 
the neutron depolarization. However, in the case of too small concentration of clusters and/or the cluster size, the 
depolarization can be also small, demanding high precision of the depolarization measurements. 

The neutron depolarization was measured at the small-angle polarized-neutron facility “Vector” (the reactor WWR- 
M, Gatchina, Russia) [IM in the temperature range 50-290 K, covering the relevant temperatures. The polarized 
neutron beam with the wavelength A = 0.85 nm and the spectrum half-width AA/A = 0.25 passed through the 
cylinder sample 3.8 mm in diameter and 13.8 mm high, that was cut from the same parent single crystal, from which 
the sample for the ac measurements was cut off. A diaphragm at the sample was 2x12 mm^. The guiding magnetic 
field 5 Oe at the sample had a minor effect on the clusters state, if at all. To ensure the absence of temperature 
hysteresis in the T-range under study, the measurements were performed both with decreasing and increasing the 


4 


temperature. The measurement time 5 h at each temperature point was sufficient for the accuracy needed. 

The linear ac MS was measured, as well, in zero dc field with the ac field of the frequency 95 KHz and the amplitude 
1 Oe in the wide temperature range 80 K < T < 180 K, above the freezing temperature. 


III. DETAILS OF DATA TREATMENT 


The experimental data treatment was performed at the PNPI computer cluster. 

The M 2 measurement data were fitted with the model function obtained from the stationary solution of FPE 
and [ 2 ^, 

dW B 

2rAr—=-u- (VH X VIT) + V(VH + /3ITVy) , (1) 

at a 

derived from the underlying Landau-Lifshitz-Gilbert stochastic equation for the magnetization vector of a single¬ 
domain particle. In Eq. W is the nonequilibrium probability-density function for directions u = Mg/Mg of the 
magnetization Mg, V = d/du is the gradient operator, t is time, V is the magnetic potential, and /3 = v/ksT, where 
V is the particle volume, and ks and T are the Boltzmann constant and the temperature, respectively. The diffusion 
relaxation time was taken in Landau-Lifshitz form as r^v = To/a, where a is the damping factor and To = 
with 7 being the gyromagnetic ratio. 

The magnetic potential V is assumed to be uniaxial [20l - l^ . viz.. 


BV = asin^ B - - ^^u-—coswt. ( 2 ) 

H h 

The first term in Eq. m is the anisotropy energy with a = B^a where Ka is the anisotropy constant and B is the 
angle between the magnetization vector and the anisotropy axis. The second and the third terms are the Zeeman 
energies of the particle magnetic moment in the bias dc magnetic field H and the ac driving magnetic field h of the 
frequency w, respectively, with = B^sH, = B^sh. No restriction on the sign of a is implied. 

In the framework of the formalism, the time scaling factor rjv in Eq. m denotes the characteristic time of diffusion 
in the absence of potential (free-diffusion time). In the case cr > 0, it is a pre-exponential factor in the conventional 
Neel relaxation time derived for a double-well potential with a barrier. 

The first and the second terms in the right-hand side of Eq. © are responsible for precession and thermal relaxation, 
respectively. When the magnetic field is parallel to the anisotropy axis, the problem simplifies considerably [s^. In 
the present study, however, this was not the case. First, by technological reasons, the anisotropy axis was oriented 
by the angle B ~ 41.4° relative to the sample plane and, hence, to the applied in-plane magnetic fields. Second, the 
point Laue patterns indicated some twinning in the single crystal used, usually observed in pseudocubic perovskites 
(84| . necessitating averaging over the twins orientations. And third, the axial symmetry, if conserved, degenerates 
precession, eliminating the precession term in Eq. o, what strongly diminishes the informative content of the data. 
An analytical solution of Eq. o for the common case of oblique magnetic field is absent, and one needs to solve 
the problem numerically. However, by expanding the solution in series by spherical harmonics, it is still possible to 
reduce the problem to an infinite system of equations, which, in turn, can be expressed as a continued-fraction matrix 
relation (see appendix). 

An experimental geometry with the parallel orientation of dc- and linearly polarized ac magnetic fields somewhat 
simplifies the data treatment, with no damage for completeness of the information content. 

The solution accuracy depends, mainly, on the number of retained equations in the system (the number of iterates 
in solving the continued-fraction equation) and on the number of retained harmonics in the Fourier expansion. The 
former number no = 8 and the latter ko = ±4 appeared to be sufficient for the whole data array. 

According to common symmetry requirements, the second harmonic must be antisymmetric on the dc magnetic 
field. However, due to the finite cycling frequency of H, hysteretic behavior may occur, violating the if-antisymmetry 
for a single, direct or reverse, iJ-scan, while the antisymmetry still conserves for the whole hysteresis loop. As the 
present stationary solution, with the steady field H, does not account for the hysteresis, an average between the direct 
and reverse scans was taken for each temperature, which, as expected, turned out to be antisymmetric within the 
experimental error. This average was additionally antisymmetrized relative to H = 0. Such elimination of a small 
hysteresis was believed not to distort the resultant parameters considerably. 

The real and imaginary parts were simultaneously fitted with the following parameters: the saturation magnetization 
of the cluster ensemble M, the anisotropy energy of a cluster Sa = VcKa, where Vc is the mean cluster volume, the 
mean magnetic moment of clusters nic = VcMg, the dispersion of the log-normal distribution Dy over cluster volumes, 
the damping factor a and two parameters describing the real and imaginary parts of the signal from the matrix 
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assumed to be linear functions of the dc magnetic field. The matrix response, when observed, includes the signals 
coming from magnetic inhomogeneities other than SPM clusters, e. g., from spin polarons. Due to their small size, 
they exhibit only slight non-linearity in weak fields and make up a minor contribution to the total response in the 
measured region H < 300 Oe, despite that these entities may occupy a large volume of the sample. 

Some additional magnetic characteristics of the cluster system can be calculated from the fit parameters, as well. 
These are the concentration of clusters N = M/rric, the mean intercluster distance (r) oc the incluster 

anisotropy field Ha = Eaj'fnci the diffusion relaxation time of the cluster magnetic moment tjv = rric/^^OikBT, and 
the characteristic dipolar energy associated with the cluster subsystem £d = 47rm^/V, where V oc N~^ is the mean 
volume per cluster evaluated from the concentration of clusters with the known perovskite-cube volume Vo ■ 

With the neutron depolarization data additionally involved, some more quantities characterizing the SPM system 
can be obtained. The polarization of the passed neutron beam can be expressed, in notations of Ref. [35l| . as 


where Po is the initial polarization directed along the beam, 7 „ is the neutron gyromagnetic ratio, Vn is the neutron 
speed, so that ^njVn ~ 46.3A nm-Oe“^, B = 47r(/r)/t)o is the magnetic induction in FM clusters with (^) being the 
mean magnetic moment per formula unit, 7?. is the mean radius of clusters, C is the weight (or volume) fraction of 
clusters in the sample and L is the sample thickness along the neutron beam. The value L = 3.58 mm was used as 
an average over the diaphragm width 2 mm. 

Equation makes it possible to separate out the incluster magnetic moment (/x), the mean cluster volume Vc-, and 
the fraction of clusters C. By definition, the mean cluster magnetic moment reads: 

TOc = — (a*) , (4) 

Vo 


noting that Vcjvo is the mean number of formula units per cluster. The saturation magnetization M, a normalization 
factor for the M 2 response function in the units emu-g“^, can be presented as 

where Na and Aimoi are the Avogadro number and the molecular weight, respectively. From Eq. one obtains: 


CVoilJif 


\QL^\^TT^ri) \P ) 


exp 



( 6 ) 


Remind that Dy is a dispersion of the volume distribution of the cluster ensemble. The exponential in Eq. is a 
correction factor to transfer from the cluster volume with the radius TZ (Eq. ([3])), vu = A-kP? ji, to the mean cluster 
volume Vc to be obtained. This radius can be expressed as an average over the volume distribution in the form 
P (X J dv-^f{v)v/v where n is a median of the assumed log-normal distribution f{v). The term v/v occurs since 
f{v)dv is taken as a fraction of the total volume occupied by the clusters with volumes in the interval (u,u + du). 
With the same note, Vc = J dvvf(v)v/v just equals to the median v. As a result, one obtains Vc = wk exp(—3Dt,/2). 
In this study, the exponential factor yields a small correction 4%-ll% to the cluster size. 

The quantities C, Vc, and (/x) can be resolved from the system of Eqs. (jH), ([S|) and ([5]). Eurther on, the mean 
cluster size will be referred to as a diameter relating to the mean cluster volume Vc = 

Thus, in this study, the neutron depolarization analysis adds to more complete quantitative characterization of the 
SPM system. 

A constant value of the perovskite-cube volume Vo = 0.0558 nm^ was used for evaluation of the parameters in the 
relevant interval 97 K < T < 140 K neglecting its relative variation of the order 1 • 10“^ due to thermal expansion 

M- 


IV. RESULTS AND DISCUSSION 

The real [Figs. l(a)-l(d)] and imaginary [Figs. l(e)-l(h)] parts of the M 2 signal are presented as functions of the 
scanned magnetic field for four characteristic temperatures. 

In Figs. 1(a) and 1(e), the signal for T = 156 K is presented as a typical example of the response in the wide range 
147 K < T < 230 K. The regular part of the signal, noticeably contaminated with the apparatus noise, is suggested 
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FIG. 1. (Color online) The second harmonic of ac response as a function of dc magnetic field for four characteristic temperatures: 
the real part (a)-(d) and the imaginary part (e)-(h). Black solid and blue open circles denote direct and reverse scans, 
respectively. The red curves on the plots (b), (c), (d), (f), (g), and (h) are best fits (see text). 
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FIG. 2. Normalized polarization of the neutron beam passed through the sample vs temperature. The solid curve is a polynomial 
fit. (Insets) “Coercive field”estimated from the real parts of M 2 (a) and resistance of the sample (b) vs temperature. 


to originate from small-size magnetic inhomogeneities, which are, most likely, spin polarons and a small portion, 
C < 10“"^, of fine FM clusters located at spatial chemical inhomogeneities formed by oxygen nonstoichiometry and 
local enrichment with Sr, favorable for the clusters nucleation. A great deal of the response still extends far beyond 
the measured iJ-field range. The marked hysteresis is due to pinning of the clusters magnetic moments. At elevated 
temperatures, the signal gradually loosens and eventually fades in the apparatus noise somewhere near T* r; 230 K 
[l|. 

At lower temperatures T < 147 K, the character of the M 2 response, exemplified in Figs. 1(b) and 1(f) for T = 141 
K, drastically changes. In a temperature interval of only 15 K, both the real and imaginary parts grow rapidly by 
an order of magnitude. The hysteretic behavior becomes much less pronounced. Upon cooling down to 140 K, the 
“coercive field”i7c2 defined by the condition ReM 2 {Hc 2 ) = 0 [Fig. 2, inset (a)] strongly falls. Such a form of the 
signal is typical for isolated SPM clusters. The extrema positions shift noticeably towards lower fields indicating the 
clusters growth, as evidenced by the calculations below. 

Starting from 130 K, the rate of the signal growth from SPM clusters strongly increases. At this stage, illustrated 
by Figs. 1(c) and 1(g) for T = 120 K, the signal becomes redistributed in favor of low fields exhibiting extrema at 
H < 50 Oe in the real part and H < 30 Oe in the imaginary part. The cluster signal dominates absolutely over 
the matrix response, the hysteresis being almost absent. This tendency persists down to the lowest temperatures 
measured. 

Below 120 K [Fig. 1(d) and 1(h) for T = 97 K], the signal continues to increase moderately with almost no change 
of its shape. This is the stage where the mean size and magnetic moment of clusters start to stabilize. 

The parameters characterizing the cluster system at temperatures 97 K < T < 140 K were evaluated from the best 
fits of M 2 {H) responses and the neutron depolarization data, as explained above. The latter are displayed on Fig. 2, 
where the polarization of the beam passed through the sample normalized by the polarization of the incident beam is 
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FIG. 3. The total saturation magnetization of the clusters ensemble (solid circles) and the characteristic dipolar energy 
associated with the cluster subsystem (open squares) (a); the mean magnetic moment rric and the magnetic moment rrim 
corresponding to a maximum position of the volume distribution (b); the concentration of clusters (solid circles) and the mean 
intercluster distance (open squares) (c); the incluster anisotropy field (solid circles) and the mean cluster anisotropy energy 
(open squares) (d), and the width of the log-normal distribution function (e) as functions of temperature. The solid curves in 
(a)-(d) are guides for the eye; the line in (e) is the linear approximation \/D^ ~ —0.349 -I- 0.0127r. 


presented as a function of temperature. The error bars are close to the size of points. The solid curve is a polynomial 
fit, to interpolate between the points. Upon lowering the temperature, depolarization of the passed beam increases, 
due to growth of the FM-cluster ensemble. At T/ « 65 K, the polarization ceases to fall and stabilizes at a constant 
level, indicating the end of FM-cluster evolution and freezing the FM-cluster dynamics on the time scale at least 10"* 
s. A value of the freezing temperature T/ is compatible with NMR data [ij and off-site ac MS measurements [s^. 
A weak inflection between 130 and 150 K correlates with the appreciable growth of FM clusters below 140 K. The 
descending temperature dependency of resistivity [Fig. 2, inset (b)] evidences a dielectric character of the sample, 
typical for the doping x < Xc- 


A. Magnetic and geometrical characteristics 

In Fig. 3, some quantities obtained solely from nonlinear magnetic response are presented as functions of temper¬ 
ature. These are the saturation magnetization of the cluster ensemble M and the characteristic dipolar energy Sd 
(a), the mean cluster magnetic moment rric and the cluster moment corresponding to the maximum position of 
the volume log-normal distribution (b), the concentration of clusters N and the mean intercluster distance (r) (c), 
the anisotropy field Ha and the mean cluster anisotropy energy (d), and the width of the log-normal distribution 
function (e). 
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Temperature dependencies of the parameters obtained by means also of the neutron polarization data are presented 
in Fig. 4, namely, the mean cluster diameter and the diameter Dm corresponding to a maximum of the volume 
distribution (a), the volume fraction of clusters C (b), and the mean magnetic moment per formula unit in the clusters 
{^i) (c). 

These dependencies specify the tendencies exemplified in Figs. l(b)-l(d) and l(f)-l(h) and visualize two stages of 
the SPM cluster MEPS alternating each other upon cooling, namely, (i) growth of FM clusters bound to chemical 
inhomogeneities in the temperature range down to 130 K and (ii) homogeneous nucleation of clusters below this 
temperature. The latter stage, in turn, proceeds via two regimes, viz., progressive increase of the mean cluster size 
down to 115 K and its stabilization below this temperature, the concentration of clusters still increasing. 

Now, these processes will be considered in more detail. At the first stage, growth of FM clusters bound to chemical 
inhomogeneities occurs gradually upon cooling [Figs. 3(b) and 4(a)]. At T r; 140 K, the clusters become appreciably 
larger than the chemical inhomogeneities. The pinning of the cluster magnetic moment strongly weakens resulting in 
reduction of magnetic hysteresis [Figs. 1(b) and 1(f)] evidenced by substantial lowering the M 2 “coercive field” down 
to Hc 2 ~ 2 Oe [Fig. 2, inset (a)]. Thus, Tg = 140 K may be accepted as the onset of SPM behavior in the cobaltite 
under study. At this temperature, the cluster magnetic moment exceeds lO^/is (150 formula units), large enough for 
the SPM formalism to become applicable. However, the clusters themselves still remain bound down to 130 K. The 
concentration of clusters and the mean intercluster distance [Fig. 3 (c)] at these temperatures vary only a little, if at 
all. The former may be considered as an estimate for the concentration of chemical inhomogeneities responsible for 
the initial cluster formation, 8.5 • 10^^ cm“^, with the mean distance (r) ~ 100 nm between them. Slow increase of 
the saturation magnetization [Fig. 3(a)] and cluster fraction [Fig. 4(b)] above 130 K occurs only due to the growth 
of the cluster size. 

At this stage, the mean magnetic moment per formula unit [Fig. 4 (c)] obeys an exponential law (/i) = 
/ie// exp(—A/T) with A = 140 K and p-e// = 2.15^b (dashed curve). Such kind of behavior, found in pure 
LaCoOs by electron spin resonance for not too high temperatures [s^j, is conventionally interpreted as thermal 
excitation to the intermediate-spin state or, alternatively, a mixture of low-spin- and high-spin states, valid for not 
too high temperatures. Thus, above 130 K, {y) of the doped cobaltite under study also follows the common tendency 
of thermal excitation, predominantly, of Co^“'" ions, implying that the incluster hole concentration only a little exceeds 
the concentration of holes in the matrix, x ~ 0.15. In this temperature region, the concentration of holes in the 
clusters is independent of temperature and determined, mainly, by the distribution of chemical inhomogeneities. 

Below T r; 130 K, the cluster fraction [Fig. 4(b)] and the saturation magnetization [Fig. 3(a)] start to rise steeply 
signifying intense growth of the sample volume involved into the SPM component. The concentration of clusters 
increases [Fig. 3(c)] manifesting entering the stage of free condensation (or homogeneous nucleation). The clusters 
are, mainly, no longer bound to chemical inhomogeneities and formation of new free clusters is activated over the 
whole sample volume resulting in increase of N by five times upon cooling to 97 K. This transition is most distinctly 
visualized via a cusp in the temperature behavior of the anisotropy field [Fig. 3(d)]. The “negative”sign of Ha 
is consistent with the “easy-plane” character of magnetocrystalline anisotropy found by neutron diffraction in this 
compound [H and . 

At the stage of free condensation, the mean magnetic moment (/r) deviates from the exponential law [Fig. 4(c)]. 
This peculiar nonmonotonous behavior of the Co magnetic moment may be explained by variation of the concentration 
of holes in clusters with temperature. At first, the concentration even slightly lowers, whereas below 115 K, holes 
start to intensely fill the clusters importing additional uncompensated electron spin s = 1/2, thus increasing (/r). Due 
to the small volume fraction of clusters, C < 10“^, this process does not cause noticeable hole depletion of the matrix. 

Two different regimes can be distinguished at the stage of free condensation. Upon cooling down to T r: 115 K, the 
magnetic moment and cluster diameter strongly increase [Figs. 3(b) and 4(a), respectively], while (/i) even slightly 
decreases [Fig. 4(c)]. Below 115 K, rric and the cluster diameters Dc and Dm cease to grow reaching the large values 
rric ~ (12 — 13) • 10^/is and Dc ~ 12.5 nm, whereas (/r) rapidly rises. The position of the distribution maximum 
stabilizes at the value Dm ~ 8 nm (compare to the long-range FM correlation length > 7 nm evaluated by Phelan 
and collaborators 0 in the polarized-neutron study of Lao.ssSro.isCoOs below 100 K). At the same time, the mean 
cluster diameter Dc starts to decrease. The latter occurs due to gradual narrowing of the (asymmetrical) volume 
distribution of the cluster ensemble from both sides towards the maximum position. From M 2 data treatment, the 
distribution width varies linearly with temperature as y^Dy r: —0.349 -I- 0.0127T [Fig. 3(e)]. Upon cooling below 115 
K, the distribution “sharpens”around the stabilized Dm- Thus, absorption of holes by clusters and increase of the 
clusters concentration are the main factors promoting the steady rise of the saturation magnetization [Fig. 3(a)] in 
the latter regime. 

Charging the clusters due to absorption of holes increases the Coulomb energy of a cluster opposing its growth. 
The alternation of the regimes at 115 K is a result of the competition between these two tendencies. 

Note also that Coulomb repulsion between the clusters may affect their space distribution especially at the lowest 
temperatures measured where the clusters are essentially enriched with holes. 
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FIG. 5. The damping factor (a) and the mean diffusion relaxation time of the cluster magnetic moment (b) as functions of 
temperature. The solid curves are guides for the eye. 


However, particular reasons for the change of the MBPS character at 130 K as well as for alternation of the regimes 
at 115 K remain an open question. 

By the way, poor convergence of the fit procedure at the high-temperature margin resulting from growing correlations 
between the parameters yields large error bars achieving 100% for some parameters at 141 K [Figs. 3(c) and 3(d)]. 
The small cluster size implies a part of the system to be somewhat beyond the SPM formalism applicability limits and 
leads to serious convergence problems. Hence, quantitative results at temperatures close to 140 K should be referred 
to with care. 
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B. Superparamagnetic dynamics 

When studying an SPM system, a conventional problem arises on the effect of interparticle interactions on SPM 
dynamics [s^. For a diluted SPM system with the insulating matrix (the present case), the most relevant is the 
long-range dipolar coupling between the particles. Its role can be estimated comparing the characteristic dipolar 
energy Sd [Fig. 3(a)] to the mean anisotropy energy of a cluster Sa [Fig. 3(d)]. The modulus of the latter well exceeds 
Ed except only for the lowest measured temperatures, thus, justifying disregard of dipolar forces in the calculations. 
At the margin T = 97 K, however, jsaj is only three times greater than Ed- Below 97 K, the cluster anisotropy energy 
tends to decrease (by modulus), while the dipolar energy increases. The ratio Ed/\Ea\ grows upon cooling signifying the 
dipolar forces coming gradually into play. Provided that this trend persists at lower temperatures, the two quantities 
meet each other somewhere above the freezing temperature Tf. In this case, the FM-cluster subsystem enters the 
regime of dipolar dynamics favoring intercluster long-range magnetic correlations and inhibiting correlations with the 
anisotropy axes directions. The SPM dynamics acquires a collective character resulting, particularly, in increasing 
the characteristic relaxation time and the blocking temperature |40l - l4^ . Eventually, the cluster subsystem freezes in 
the dipolar-cluster-glass (or super spin-glass [dlj) state. 

The damping factor a and the relaxation time ryr (Fig. 5) are also deduced only from M 2 measurements. The former 
is a purely dynamical quantity, hardly measurable with other techniques. This parameter is responsible for relaxation 
of the FM-cluster magnetic moment due to its interaction with the environment. This might be incluster lattice 
and magnetic excitations, magnetic inhomogeneities, etc. To our knowledge, there is no microscopic theory of SPM 
relaxation. In the formalism, a was assumed to be explicitly dependent neither on magnetic field, nor on any cluster- 
system parameter. Its smooth rise with temperature is characteristic of a quantity depending on thermal excitations. 
The steeper increase with temperature above 130 K may be due to the growing effect of magnetic inhomogeneities as 
the clusters become bound and their magnetic moments pinned to chemical inhomogeneities as their size decreases. 
The steeper decrease of a below 115 K is suggestive of an onset of the freezing tendency enhanced (or provoked) by the 
intercluster dipolar forces. The effect of dipolar forces on the quantity characterizing free diffusion should not confuse, 
as the latter was defined with respect to the magnetic potential ([2]) not containing intercluster interactions. Explicit 
inclusion of the dipolar-coupling term into Eq. would have eliminated this effect from a and r^r by transferring it 
to the potential-dependent part, such as an exponential in the Neel relaxation time in its conventional definition in 
the case ct > 0 . 

A magnitude of the damping factor itself yields essential information concerning the influence of precession on 
relaxation of the cluster magnetic moment. The first and the second terms in the right-hand side of Eq. o describe 
precession and thermal, diffusion-type, relaxation of the magnetic moment, respectively. In the case of axial symmetry, 
when the magnetic field is parallel to the anisotropy axis, as well as in the limit a ^ 1 (overdamped case), the first 
(precession) term vanishes in favor of diffusion. The measured values a = 0.3 — 0.6 point out that relaxation is 
considerably modified by precession. Both the terms of FPE turn to be comparable generating interplay of precession 
and thermal diffusion. Upon cooling, the role of precession increases. This phenomenon was examined in simulations 
by direct numerical solution of the Landau-Lifshitz-Gilbert equation (Ref. [331 ] and references therein). 

The relaxation time rv [Fig. 5(b)] experiences intense increase upon cooling following the combined effect of TOc, 
the damping factor, and the temperature, with the main contribution from the growing cluster magnetic moment. 
The values tm of the order 10“^° s are quite typical of SPM dynamics. The steeper growth below 115 K is suggestive 
of the beginning of the slowing-down tendency leading, eventually, to freezing the cluster system. 


C. Linear and third-order susceptibilities 

With known parameters of the SPM system, it is possible to retrieve the linear magnetic response of FM clusters 
and to estimate its role in magnetism of the cobaltite under study. In Fig. 6, temperature dependencies of the total 
magnetic susceptibility measured at 95 KHz and zero dc magnetic field and the linear SPM response recovered to these 
experimental conditions from M 2 measurements are presented. The recovered SPM susceptibility arises T k, 140 
K correlating with the onset of FM clusters growth [Figs. 3(b) and 4(a)] and SPM behavior. Upon cooling, its 
relative contribution increases and becomes dominant below 120 K, in the region of intense growth of the saturation 
magnetization M [Fig. 3(a)]. Recall, that for all temperatures measured, the clusters occupy not more than 2 • 10“^ 
of the sample volume [Fig. 4(b)], indicating that the large cluster magnetic moment is the main factor highlighting 
the SPM contribution to magnetic susceptibility. 

In the inset to Fig. 6, the third harmonic is also displayed as retrieved for zero dc field and the ac field with the 
frequency 30 kHz and the amplitude 1 Oe. Its temperature dependency with the broad maximum at T Ri 115 K 
stands alongside with the behavior of the low-frequency third harmonics measured by Caciuffo and collaborators on 


13 



80 100 120 140 160 


T [K] 


FIG. 6. The measured total linear magnetic susceptibility (solid circles and line) and the clusters linear susceptibility as 
recovered from the second harmonic data for the excitation frequency 95 kHz (open circles and dashed line). (Inset) Third 
harmonic of the response vs temperature recovered for the excitation frequency 30 kHz. The curves are guides for the eye. 


the powder La.rSri_,T,CoOs for a set of doping values in the “metallic”region where FM clusters distinctly exhibited 
themselves in a pure form [iGj. These authors observed qualitatively similar dependencies with the maxima at T Ri 175 
and 150 K for a: = 0.25 and 0.20, respectively, correlating with the maxima in respective linear susceptibilities. For 
X = 0.15, the linear response exhibited a wide hump in the interval 150 K < T < 200 K associated with FM clusters 
(unfortunately, the third harmonic was not presented for this doping). The latter is absent in Fig. 6 due to much 
lesser concentration of clusters in the single crystal. 


V. SPM DYNAMICS ANSATZ: GILBERT vs LANDAU-LIFSHITZ 


High accuracy and representativity of the experimental data enables to put a question concerning the main alter¬ 
native underlying the formalism used. 

Historically, the Gilbert dynamical equation describing evolution of the classical magnetic moment m of a large 
FM cluster was introduced phenomenologically as a starting point (unpublished work, mentioned in Ref. [i^) [s^: 


dm 

-= am X 

dt 


Be// - 


a dm 
ym dt 


( 7 ) 


Be// = - 


dV 

i9m 


where the effective field is given by 
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TABLE I. Normalized chi-squares of the best fits, the damping factor a, the relaxation time rjv, and the anisotropy field Ha 
evaluated for a set of temperatures for the cases of Landau-Lifshitz- and Gilbert schemes. 


T [K] 

\norm 

Landau-Lifshitz: tn = 
a Tn [10^° s] 

Ha [Oe] 

\norm 

Gilbert: 

a 

tn = To{a + a ^) 

tn [10“ s] 

Ha [Oe] 

97 

2.88 

0.33 

7.2 

-17 

2.33 

0.36 

6.5 

-21 

100 

2.96 

0.37 

6.2 

-20 

2.25 

0.41 

5.5 

-26 

104 

2.07 

0.40 

5.4 

-23 

1.48 

0.44 

4.9 

-33 

109 

1.94 

0.40 

5.6 

-26 

1.53 

0.42 

5.6 

-40 


with V from Eq. Equation © contains time derivatives of the magnetic moment in both the right- and left-hand 
sides. This inconvenience was eliminated by expressing the derivative explicitly, resulting in the Landau-Lifshitz 
equation [H and 13 • 


~ “ r -D 1 

-= ym X Beff — 7—m x m x Beffl 

dt ' m 


( 8 ) 


with the renormalized gyromagnetic ratio 7 = 7/(1 + a^). The famous Landau-Lifshitz relaxation (damping) term 
proportional to —m x [m x Beff] drives m to the direction of Be//, while a measures the magnitude of the relaxation 
term relative to the gyromagnetic term. 

In this transformation, no restrictions or additional assumptions were involved; thus, both the equations, © and 
are of one and the same level of generality. Moreover, both of them lead, formally, to one and the same FPE; 
hence, each of the equations may be considered as generating one for the FPE. At first sight, the choice of the 
commencing equation seems to be a scholastic question. However, due to renormalization of 7, the characteristic 
diffusion time tn entering Eq. © is different in these two cases, depending on whether Gilbert- or Landau-Lifshitz 
equation is chosen as commencing one. The diffusion time is tjv = To {a + C(~^) in the former and rjv = Toa~^ in the 
latter case. A question arises, in which scheme, Gilbert or Landau-Lifshitz, 7 has to be taken as bare one. Generally, 
this difference may result in dissimilar parameter meanings evaluated from best-fit calculations. 

One might try to distinguish between the two schemes and to make a choice in favor of one of them by comparing 
chi-squares of the experimental-data best fits. As mentioned above, in the overdamped case (large a), the precession 
term in Eq. © eliminates and the solution depends on the damping factor only implicitly via tn. Thus, one can 
hardly expect any meaningful difference in the chi-squares. In the opposite case, a <C I, there is also no difference, 
as Tjy is the same for both the schemes. The presently obtained values of a [Fig. 5(a)] lie just in the intermediate 
region for which the search of distinction in the best-fit quality makes sense. 

Appropriate calculations were performed for a number of points in the lower-temperature region where the SPM 
system is well formed, with the steady fit convergence. The results, including the normalized chi-squares and the 
most sensitive parameters among those determined solely from the nonlinear measurements, viz., the damping factor, 
the relaxation time, and the anisotropy field, are presented in the table. It is seen that the Gilbert scheme yields 
systematically smaller values of Xnorm than these for the opposite, Landau-Lifshitz, one. Thus, the Gilbert approach 
seems to be more realistic, at least, for the object under study. The two parameters shown in the table differ essentially 
for the two cases, while the differences between other parameters (not presented) are of the order 1%. At any rate, 
calculations performed within either of the alternative schemes would result in qualitatively similar temperature 
dependencies. 

To establish the extent of universality of this finding, a study of the wide set of more conventional SPM systems is 
needed. 


VI. CONCLUSION 

The study of Lao.ssSro.isGoOa in the wide temperature range from 97 to 230 K was performed with linear and 
nonlinear ac magnetic response as well as neutron depolarization techniques. The system of large ferromagnetic 
clusters was found to emerge as a particular type of magnetoelectronic phase separation, besides well-established spin 
polarons. The experimental data were treated with a thorough formalism based on the solution of the Fokker-Planck 
equation to obtain magnetic, geometrical, and dynamical parameters characterizing the cluster ensemble. 

Evolution of the cluster system upon cooling was found to occur via two processes, namely, (i) gradual growth 
of clusters bound to local chemical inhomogeneities in the wide temperature range down to 130 K with reaching 







15 


the superparamagnetic regime at 140 K accompanied by failure of pinning the cluster magnetic moment and (ii) free 
condensation (or homogeneous nucleation) below 130 K consisting in intense formation and growth of clusters over the 
whole sample volume accompanied by variation of the hole density in the clusters. The latter process proceeds in two 
regimes alternating each other at 115 K, viz., growth of the cluster size with the incluster hole content only slightly 
exceeding the nominal doping x = 0.15 above this temperature and stabilization of the cluster size accompanied by 
intense increase of the magnetic moment per unit cell due to enrichment of clusters with holes below 115 K. 

At the onset of superparamagnetic behavior, ferromagnetic clusters start to contribute noticeably to the total linear 
magnetic susceptibility and, below 120 K, their contribution becomes dominant, at least down to 97 K. 

The growing effect of intercluster dipolar coupling upon cooling may result in a dipolar character of the frozen state 
below 65 K. 

The observed ferromagnetic clusters are suggested to be (i) a precursor of the ferromagnetic state at higher doping 
levels and (ii) an onset of large ferromagnetic clusters revealed by neutron scattering below 100 K Q. Thus, the 
phase separation and its evolution in a single crystal at temperatures higher than 100 K was shown to have a more 
complicated character than was supposed before. 

Relaxation dynamics of the cluster magnetic moment was found to be, essentially, the interplay of precession and 
thermal diffusion. Below 115 K, the magnetic dynamics exhibits a tendency to freezing favored by intercluster dipolar 
coupling. 

High representativity and accuracy of the ac measurements performed made it possible to elucidate the ground 
alternative underlying superparamagnetic dynamics and to make a choice, though not robust enough, in favor of 
Gilbert approach. 

The method employed was demonstrated to be an efficient means for studying the systems exhibiting superparam¬ 
agnetic behavior. 
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Appendix 


The formalism used is presented in Refs, 
series over spherical harmonics. 


2(1 l2ll . and |2J . By expanding the distribution function in Eq. o in the 


00 I 

Wif) = EE C/m (^)Em j 

m— — l 

the problem is reduced to the recurrence relation and |2l| , 

Sn = —[Qn + Q^S„_|_iQ„_|_i] ^ , 


(A.l) 


(A.2) 


relative to the matrix continued fraction S„. The tridiagonal supermatrices Q^, given by 

[Qn]/,m — —l,mPn T T : 

where q„(ma;) = —irriTNUjl + q„, with I being the identity matrix. The supermatrices p„, p^, q„, and q(^ are 
presented as follows poj| : 


P- Oj 


Pn = 


f ^2n 

\ b2„_l 


q 


+ 

n 


Z2n 

0 


d2„ \ 
a2n-l J 

Y2„ \ 

Z2n-1/ 


Pn 


fo b2n\ 

lo 0 ) 


J qn 

J qn 


f V2n 0 \ 

{W2n-1 V 2 U- 1 ) ’ 

} X2„ W2„ ^ 

\Y2n-l X2n-1 J 


Dimensions of the supermatrices p„ (qn), Pn (qli), p„ (q„) are 8n x 8n, 8n x 8{n + 1) and 8n x 8(n — 1), 
respectively. Their matrix elements are, in turn, the matrices with the elements given by 
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(^/)n,m — —l,m^/ “1“ ^n,7nO>l, — l-\-7n — l _/_^77j_2 : 

(t>/)n,m — ^71,771^1 —I + 5n+l ,m —/+m—2 H“ ^n+2,m^/^_/_^yjT,_3 ; 

(d/)7T,^77T, — (^n—2,m^/ _^_|_772-(_i ~l~ ^n——/+m “t“ _/_|_772_i ; 

— ^n— —iJ^Yn ^n,m^l, — l-\-m — l ^n-\-l,m^i _/_|_7^_2 : 

(Y/)yj 772 — _/_|_772_i “1“ ^n-\-l,myi.—l-\-m —2 H“ ^n+ 2 ,my/ _/_^772_3 j 

{y^l)n^m — ^n—2,m'^^/^_/_^772+l /+m H“ ^n,m'^/^_/_|_ 772 _i ; 

(^/)n,m — ^n-\-2^m^l, — l-\-m— 3 j ('^[}n,m — ^Ti—2,m‘l^l, — l+m+l ■ 

The last five supermatrices are drawn explicitly in Ref. [^ . 

Using the direction cosines for the vectors H and h, 


7i = sin Ip cos (p , 72 = sin ip sin cp , 73 = cos ip 


and 


7 ' = sin ip' cos (p' , 72 = sin ip' sin (p' , 73 = cos ip' , 

respectively, the matrix elements in Eqs. (IA.3I) can be written in the form [lO and [24 

,"1473 


(^n.m — 


4 a 


+ .4(71-42) / 7 —^-TUTT-T 
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8a 
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4 Y (2n+l)(2n-l) ’ “ 8 y (2n+l)(2n-l) 


(T[n(?T. + 1 ) — 3to^] n(n+l) 

TT7TZ ^ 


(2n- l)(2n + 3) 


2a 


+ .^//(7i-42) / 7 —^-rUTT-T 

xZ^ = —i -^-V (" + w + lj(" — fn) 


4 a 


(inl^n ^ .ctto^ / (n + 1)2 — TO^ _ ^ (71 “ 42 )"- /(n + to + l)(n + to + 2 ) 

~ *44y y ( 2 n + l)( 2 n + 3) ’ “ 4 V ( 2 n + l)( 2 n + 3) 


, n + 1 CTTO 

W„,m = -*- 

' 2 a 


'n? — m? I (71 — 42 )(" + 1 ) {n — m){n — m — 1) 

44-1 ’ ^ 4 


4n2 — 1 


^n.m — 


CTO / [(n + 2)2 — TO2][(n + 1)2 — to2] 


2n + 3 y 
with the additional relations, 
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The stationary ac response can be calculated from the continued-fraction solution Si of Eq. (|A.2p . viz., 
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where each element of Ci is a column vector, 


and 
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(A.4) 


Elements of the column vector cf(a;) are Fourier components of the corresponding coefficients cim in the expansion 

O. rn r-i 

Following Refs. and[^, we define the ac response as magnetization Mh{t) in the direction of the driving field 
h, with the Fourier transform, 

00 
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with o(u;), _i(w), and ^(w) taken from Eq. (IA.4I) . 
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